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ABSTRACT 



Context. The study of pulsation in pre-main-sequence intermediate mass stars represents an important tool to derive information 
q , on the stellar parameters and structure, as well as to test the validity of current theoretical models. The interest for this class 
$_i ' of variable stars has significantly increased during the last decade and about 30 members are presently known in the literature. 
Aims. A new observational study of the Herbig Ae star VV Ser has been performed in order to detect and accurately measure 
■ pulsation frequencies in the S Scuti range, thus enlarging the sample of known pulsators and contributing to the empirical 
^ definition of the pre-main-sequence instability strip. As it belongs to the continuous field of view of the asteroseismological 
• i-h . satellite COROT, this study also aims at characterizing the properties of VV Ser as a potential "COROT additional program" 
1 candidate. 

' Methods. CCD time series photometry in the Johnson V filter has been obtained during three consecutive years. The resulting 
light curves have been subject to detailed frequency analysis and the derived frequencies have been compared to model predic- 
tions. 

Results. Seven pulsation frequencies have been measured on the basis of the best dataset obtained in 2004, ranging from ~ 31 to 
~ 118 fj,Hz, with an accuracy of the order of 0.5 ^Hz. The comparison with an extensive set of asteroseismological models shows 
that all the observed periodicities can be reproduced if the stellar mass is close to 4Mq. Conversely, the measured frequencies 
can be associated to p modes only if the effective temperature is significantly lower than that obtained from the spectral type 
conversion. 

Conclusions. The present results seem to suggest that more accurate spectral type determination is necessary in order to dis- 
criminate the best fit model solution. In any case, the stellar mass of VV Ser is close to the upper mass limit (~ 4M©) for this 
class of pulsators. 

Key words, stars: variables: 5 Scuti - stars: oscillations - stars: pre-main sequence - stars: fundamental parameters - stars: 
individual VV Ser 



1. Introduction 

Pre-Main sequence (PMS) stars with masses larger than 
1.5 Mq are known as Herbig Ae/Be stars ( |Herbig 1960| ). 
In general, they are found within star forming regions 
and show variable emission lines (especially Ha), as 
well as strong infrared excess caused by the presence 
of circumstellar material (dust). In addition, Herbig 

Send offprint requests to: V. Ripepi 



Ae/Be stars are characterized by photometric and spec- 
troscopic variability on time scales of minutes to years, 
mainly due to photospheric activity and interaction with 
the circumstellar environment (see e.g. lOahm et al. 1 995 
IBohm fc (Mala~T995ll. 

Considerable theoretical work has been done recently 
that has advanced our understanding of PMS evolution 
(e.g. IPalla fc Stahler 19931 ID'Antona fc Mazzitelli 19941 
ISwenson et al. 1994(1 . Yet, there remain differences in the 
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models owing to alternative treatments of convection, 
opacities and the zero-point of the calculated ages. It is 
therefore desirable to find independent tools to constrain 
PMS evolutionary tracks and, in turn, the internal struc- 
ture of intermediate-mass PMS stars. 

Asteroseismology of Herbig Ae/Be stars can in prin- 
ciple allow us to test PMS models by probing their 
interiors. It is now well established that intermediate- 
mass PMS stars during their contraction towards the 
MS cross the pulsation instability strip of more evolved 
stars, giving origin to a variable class named PMS 5 
Scuti stars (Kurtz & Marang, 1995, Marconi & Palla, 
1998, ICatala 200.11 |Ripepi et al. 2006| ). The existence of 
pulsating Herbig Ae stars was originally proposed by 



Breger (1972) who discovered two candidates in the young 



open cluster NGC 2264. More than 20 years later, the 



suggestion was confirmed by Kurtz & Marang (1995) and 



Donati et al. (1997) who observed 6 Scuti-like pulsations 
in the Herbig Ae stars HR 5999 and HD 104237, respec- 
tively. 

The first theoretical investigation of the PMS instabil- 
ity strip based on non-linear convective hydrodynamical 



models was carried out by Marconi & Palla (1998) who 
calculated its topology for the first three radial mod es. A 
subsequent theoretical work by Suran et al. (2001) made 
a comparative study of the seismology of a 1.8 M PMS 
and post-MS star. Suran et al. (2001) found that the un- 
stable frequency range is roughly the same for PMS and 
post-MS stars, but that some nonradial modes are very 
sensitive to the deep internal structure. In particular, it 
is possible to discriminate between the PMS and post-MS 
phase using differences in the oscillation frequency distri- 
bution in the low frequency range (g modes). 

Since the work by Marconi & Palla (1998) many 
new PMS S Scuti candidates have been observed and 
the current number of known or suspected candidates 
amounts to about 30 stars (see e.g. IZwintz et al. 20041 
|Ripepi et al. 2006"| ). However, only a few stars have 
been studied in 



detail (e.g. 


Ripepi et al. 2003 


IRipepi et al. 2006| 


, so that the 



overall properties of this class of variables are still poorly 
determined. 

In this context, our group has started in the year 
2000 a systematic photometric monitoring program 
of Herbig Ae stars with spectral types from A to F2- 
3 (see e.g. IMarconi et al. 20001 |Ripepi et al. 2002| 
IPinheiro et al. 20031 |Ripepi et al. 2003| 

|Ripepi et al. 2006"| ). The aims were: 1) to identify 
observationally the boundaries of the instability strip for 
PMS S Scuti pulsation; 2) to study in detail through 
multisite campaigns selected objects showing multiperi- 
odicity. As part of this long-term project, we present 
here a new observational study devoted to the Herbig Ae 
star VV Ser. In addition to the motivation mentioned 
above, this star conveniently falls in the continuous 
field of view of the asteroseismological satellite COROT 
(corot.oamp.fr). Thus, VV Ser represents a potential 
"additional program" candidate ( W eiss et al. 2 004). 



Table 1. Literature values of Spectral Type, Luminosity 
Class and rotational velocities for VV Ser 



Source 



Sp. Type 


vsini 




km/s 


A2e 




A: 




Ble-B3e 




A2eV/3 




B96 1 


200 2 


B5e 






85 3 


A3IIe/? 




AOVevp 4 


229 ± 9 5 


B+sh 


142 


A2IIIe 6 




B6e 





Herbig (1960) 
|Cohen fc Kuhi (1979)1 
|Finkenzeller fc Mundt (1984)| 
|Chavarria et al. (1988)| 



|Hillenbrand et al. ("1 992) 
|Fernandez et al. (1995)| 

|Grady et al. (1996j| 
|Gray fc Corbally (199871 
|Mora et al. (200"T)]~ 
|Vieira et al. (2003)| 
|Acke fc van den Ancker(2004j1 
|Hernandez et al. (2004)| 



l=photometrically determined 

2=no reference is indicated for this value 

3=the Authors quote Hillenbrand et al. (1995, private commu- 
nication) as source for this value 
4=the Authors quote an error of five classes 
5=vsini derived using just one line (Mgll 4481 A line) 
6=the Authors quote van den Ancker (private communication) 
as source for this value 

Table 2. Stellar parameters for VV Ser found in the lit- 



erature 


V 


D 


A v 


l0g(T eff ) 


log(L) 


M 


source 


mag 


pc 


mag 


K 


Lq 


M 




12.666 


245 


6.1 


3.95 


> 1.62 




1 


11.87 


440 


3.0 


4.03 


1.8 


3.3 


2 


11.92 


440 


3.4 


4.14 


2.23 


3.8 


3 


11.63 


296 


3.4 


3.95 


1.51 


2.1 


4 


11.58 


440 




4.03 


2.03 




5 




330 


2.7 


3.95 


1.27 
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l=| Chavarria et al. (1988)|2 = |Hillenbrand et al. (1992)| 
3= |Hernandez et al. (2004)| 4=Rostopchina et al. (2001) ; 
5=Testi et al. (1998); 6= |Acke fc van den Ancker(2004)| 



In the next section, we discuss the properties of VV 
Ser; in Section 3 we describe the observations and data 
reduction; in Section 4 we deal with the frequency analy- 
sis; in Section 5 we present a comparison with theoretical 
predictions. Finally, a brief summary closes the paper. 

2. Stellar parameters of W Ser 

The variable star VV Ser was identified as a young Ae 
star in the seminal work by Herbig (1960) on emission 
line stars associated with reflection nebulae. The star is 
located in the Serpens molecular cloud (Chavarria et al. 
1988) and has been widely studied in literature, although 
its properties and position in the HR diagram are still 
rather uncertain . In the following, we present an overview 
of the current observational status of VV Ser. 

— Spectral type and luminosity class: We report 
in Table ^ (first column) all the independent mea- 
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surements of spectral type (ST) found in literature. 
The uncertainty is large with ST varying between 
Ble-B3e and A3e. The determination of the earliest 



type by Finkenzeller & Mundt (1984) is based on the 
Hel 5876 A line which was later found to be as- 
sociated with the hot regions of accretion disks (see 
Rost opchina et al. 200 1| and references therein) and 
therefore should not be considered reliable. However, 
a B5-B9 type has also been claimed by other authors 
(see Table^). On the other hand, Table^hsts six stud- 
ies where VV Ser is classified as an Ae star, with typical 
value A2. Interestingly, the only paper exclusively ded- 



icated to VV Ser is that by Chavarria et al. (1988) who 
assigns an A2e class, based on high resolution spec- 
troscopy and Stromgren photometry. In the following, 
we will rely on these results, allowing an uncertainty of 
2 subclasses to take into account the spread of spectral 
types found in the literature. 

Concerning the luminosity class of VV Ser, there are 
three determinations in the literature: class II (Gray 
& Corbally 1998), III (Acke & van den Ancker 2004), 
and V (Mora et al. 2001). Again, there is no agreement 
between different authors. 

Distance and position in the HR Diagram: 

The distance estimate of VV Ser is related to that 
of the Serpens Cloud. On the basis of photomet- 
ric and spectroscopic observations of 5 stars in 



the cloud (excluding VV Ser), Chavarria et al. (1988) 



and de Lara et al. (1991) calculated a distance of 
D=311±38 pc, in fair agreement with previous de- 



terminations of about 440 pc by Racine (1968) and 



Strom (1974) In the following, we will adopt a dis- 



tance of 300-400 pc. 

As for the position in the HR diagram, Table lists 
several estimates of log(T ff) and log(L/L ) found in 
literature. The table highlights the existing uncertain- 
ties, mainly because of different assumptions on the 
ST, on the visual absorption and on the apparent vi- 
sual magnitude of VV Ser which varies by tenth of 
magnitude on time scales of weeks. 
Rotational velocity: Another source of uncertainty 
in the analysis of VV Ser is represented by its 
rotational velocity. Table Q reports four indepen- 
dent measurements of wsini found in the literature. 



The values reported by Hillenbrand et al. (1992) and 



Mora et al. (2001) are in fair agreement, whereas those 



given by Grady et al. (1996) and Vieira et al. (2003) 



are significantly different. More in detail, we were 
not able to find the source of the v sin i reported by 
Hille nbrand et al. (1992)| who quoted the studies of 



Finkenzeller (1985)| and |Davis, Strom fc Strom (1983) 
that however do not give a value for VV Ser. Similarly, 
the low vsini reported by Gray & Corbally (1998) 
cannot be verified because it is quoted as pri- 
vate communication. Finally, the v sin i measured by 



Mora et al. (2001) is somewhat uncertain because it 



relies on just one line, the Mgll 4481 A line 



Thus, the rotational velocity of VV Ser is still uncer- 
tain, even if there are indications that it could be larger 
than 100 km/s. This possibility will be taken into ac- 
count in the interpretation of the data (see Section 5) . 

In conclusion, the empirical estimate of the stellar 
parameters of VV Ser is still plagued by large errors. 
However, as we will discuss in Sect. 5, the pulsational 
analysis of the detected frequencies can help to reduce 
the uncertainty on luminosity, effective temperature and 
mass. 

3. Observations and data reduction 

The observations of VV Ser have been carried out in sev- 
eral runs during 2002, 2003 and 2004 (a detailed obser- 
vational log is given in electronic form in Tabled, us- 
ing the Loiano 1.5m telescope equipped with the BFOSC 
instrument (see |http://www.bo.astro.it/loiano/index.htm| 
for a detailed description of the telescope and of the in- 
struments). The CCD is a EEV 1340x1300 pixels with 
individual size of 0.58 arcsec, for a total field of view of 
13 x 13 . In total, about 55 hours of observations have 
been gathered in the Johnson V filter. 

The comparison star used was USNO-A2.0 0900- 
12930815 (18 h 28™ 33.535 s +00° 02' 25. 7l" J2000) which 
is about 7.3™ away from VV Ser, allowing to observe 
both stars in a single BFOSC frame. We note that this 
star has already b een used as comparison for VV Ser by 
Fernandez (1995) who found it stable at a level of 0.02 mag 
and with an apparent magnitude and color of V = 12.28 
mag (B — V) = 0.89 mag, respectively. 

All the data have been reduced following the usual pro- 
cedures (de-biasing, flat-fielding) by using standard IRAF 
routines. The aperture photometry has been carried out 
using routines written in the MIDAS environment. The 
typical precision of the data was of about 4, 6, and 3 mmag 
for the 2002, 2003 and 2004 dataset, respectively. On this 
basis (see also Fig.^ an d |2J it is evident that the 2002 and 
especially the 2004 datasets are the best ones, whereas the 
2003 run is of poor quality. Note that, because of the high 
sampling, the 2004 data have been smoothed by using a 
3-point boxcar filter. 

The observed light curves for the different datasets are 
reported in Fig. ^ It is evident that VV Ser displays mid- 
term variations on time scale of days, in addition to the 
pulsation on time scale of hours. In order to eliminate 
these night-to-night variations and to prepare the data 
for Fourier analysis, we decided to detrend the data to a 
common average zero value. Then, we used the period04 
package ( |Lenz &: Breger 2005| ) to remove the frequencies 
0.60 c/d from the 2002 dataset, 0.46 c/d+0.19 c/d from 
the 2003 dataset, and 0.45 c/d from the 2004 dataset, 
with the exception of the two nights with HJD=53209 
and 53210 for which a simple average was subtracted. As 
a result, we show in Fig. the light curves for the three 
datasets after the detrending procedure. 
Before proceeding with the frequency analysis of the ob- 
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Fig. 2. Left panels: light curves for the 2002 (top) and 2003 (bottom) datasets after the detrending. Note that 
SV = Vvar — Vcomp- Right panels: as before, but for the 2004 dataset. In all the panels, the solid line displays the 
fit to the data with all the significant frequencies found for each dataset, as listed in Table [3] 
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Fig. 1. From top to bottom, the first two panels show the 
2002 and 2003 dataset; the third and fourth panel show 
the 2004 dataset. Note that in the latter we have added 
0.03 mag to the original differential photometry (see Sect. 
3). The solid line shows the fit with a sin curve with fixed 
period of 2.2 d and free amplitude and phase (for each 
dataset, see Sect. 3). 

tained time series, we would like to discuss more in detail 
the mid-term light variations with period of the order of 2 
d which have been removed from the data as described be- 



fore. To investigate further this point, we fitted a sinusoid 
with the period P=2.2 d (frequency 0.45 c/d) found using 
the 2004 data (HJD=53182-53186) to the three datasets, 
leaving the amplitude and phases as free parameters (to 
this aim we added 0.03 mag to the photometry obtained 
during the nights with HJD=53209 and 53210). The re- 
sult is shown in Fig.^by a solid line: the fit is fairly good 
also for the 2002 dataset, whereas it is worse for the 2003 
dataset. In conclusion, a periodicity of the order of 2 d 
seems to be present in VV Ser. It is worth to note that 
analysing six year data Shevchenko al. (1993)| found both 
a long-term periodicity (around 1000 d) and a short-term 
periodicity (which they quote as uncertain) of 3.9 d. Even 
if there is a discrepancy between Shevchenko al. (1993) s 
results and ours, the short-term periodicity could be real. 
The physical origin of this kind of periodicity is not eas- 
ily understood and in any case is not the subject of this 
paper. However, such short-term variability in Herbig Ae 
stars is generally attributed to rotational modulation due 
to magnetic activity (see e.g. lCatala et al. 19 991. 



4. Frequency analysis 

The detrended time series shown in Fig. |21 have 
been analysed by using again the period04 package 
| |Lenz fc Breger 2005| l, based on the Fourier transform 
method. For a better interpretation of the results, we have 
first calculated the spectral window (SW) for each dataset. 
The result is shown in Fig. where from top to bottom 
we report the SW for the three time series identified in 
the previous section (see labels in the figure). The SW 
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-4 -2 2 4 
Frequency (c/d) 



Fig. 3. Spectral window in amplitude for the three data 
sets analysed in this paper (see labels in the figure). 

Table 3. Frequencies, amplitudes and confidence levels 
for the three datasets analysed in this paper. The error on 
the 2002, 2003, and 2004 data are -0.14 c/d, -0.22 c/d, 
~ 0.04 c/d, respectively. Note that for the 2002 and 2003 
datasets the frequencies are labeled in correspondence to 
the 2004 dataset (see text). 

Frequency Amplitude Confidence 





(c/d) 


(mmag) 


(%) 






2002 DATA 




h 


5.11 


6.3 


99.9 


h 


3.84 


4.6 


99.9 


k 


9.55 


2.9 


99.9 






2003 DATA 




h 


3.99 


7.6 


99.9 


k 


8.50 


2.3 


99.0 


h 


6.11 


3.1 


90.0 






2004 DATA 




h 


6.12 


6.4 


99.9 


h 


2.69 


7.8 


99.9 


k 


4.47 


3.7 


99.9 


h 


3.90 


4.8 


99.9 


k 


9.96 


2.7 


99.9 


k 


7.56 


2.1 


99.9 


h 


10.24 


1.4 


90.0 



whereas Fig. 0] (right panels) displays the 2004 one. Here, 
in each panel the peak with largest amplitude is selected 
and then removed, obtaining a new spectrum shown in 
the next panel. The last panel shows the periodogram 
after the prewithening with all the significant frequen- 
cies. The solid, dashed and dotted lines represent the 
99.9%, 99% and 90% confidence levels calculated follow- 



ing the widely used r ecipe by Breger et a l (1993) and 
Kuschnig et al. (1997) The error on the measured fre- 
quencies (apart from the ± 1 c/d alias) can be roughly es- 
timated from the FWHM of the main lobe in the spectral 
window (see I Alvarez et al. 19981 and references therein) . 
As a result, we found A/ ~ 0.14 c/d, 0.22 c/d and 0.04 
c/d for the 2002, 2003 and 2004 datasets. For the last 
dataset, due to the complex structure of the main lobe of 
the spectral window, in order to be conservative, we dou- 
bled the nominal error. It is important to emphasize that 
all the frequencies reported in the following are affected 
by the indetermination due to the 1 c/d alias. 

All the significant frequencies found for the three 
datasets are summarized in Table |3| A fit to each dataset 
with the quoted frequencies is shown in Fig. by a solid 
line. 

In order to discuss in detail the results of the frequency 
analysis, we will take as reference the frequencies obtained 
with the best dataset, that of 2004. Then, we find: 

— fx = 6.12 c/d: it is found in all datasets within the 
quoted errors. The first frequency in 2002 is clearly 
the -1 c/d alias (this periodicity has also the same am- 
plitude as in 2004), whereas in 2003 it has a signifi- 
cantly lower amplitude, owing to the poor quality of 
the dataset and/or to a possible change in amplitude 
(common occurrence among S Scuti stars) 

— fii hi f&, h- they have been found only in the 2004 
dataset. We note that fj would be exactly twice /l if 
we took the -1 c/d alias, therefore it could be a non- 
independent frequency. Moreover it is only partially 
significant (90%). 

— fi = 3.90 c/d : within the errors, it can be considered 
the same frequency as in the 2002 (/ 4 =3.84 c/d) and 
2003 (/ 4 =3.99 c/d) datasets. 

— f 5 = 9.96 c/d: it could correspond to the frequencies 
at 9.55 c/d and 8.50 c/d (-1 c/d alias) in the 2002 
and 2003 datasets. Although these three frequencies 
are not equivalent within the errors, we will assume 
that they represent the same frequency because of the 
similar amplitudes. 

In the following, we will use the set of 6 significant fre- 
quencies extracted from the 2004 dataset for comparison 
with theoretical models. 



was used as a diagnostic to distinguish between real and 
spurious frequencies. 

Each data set described in the previous Section has 
been analysed separately. Figure 0] (left panels) shows 
the Fourier transform for the 2002 and 2003 datasets, 



5. Comparison with theory 

We now attempt to use the frequency data reported above 
to carry out a preliminary comparison of the available ob- 
servables with theoretical models. To do so, we consider 
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Frequency (c/d) Frequency (c/d) Frequency (c/d) Frequency (c/d) 

Fig. 4. Left panels: Frequency analysis for the 2002 and 2003 datasets. The solid, dashed and dotted lines show the 
99.9%, 99% and 90% significant levels. In each panel, one peak (i.e. the labelled frequency) is selected and removed 
from the time series and a new spectrum is obtained. The last panel displays the periodogram after the prewithening 
with all the significant frequencies. Right panels: As before, but for the 2004 dataset. 



the observed location of the star in the theoretical HR di- 
agram and how a tentative mode identification of the ob- 
served frequencies can constrain the stellar mass. Because 
the outcome depends strongly on the assumptions, we will 
also evaluate the effects of the uncertainties on the results. 

5.1. Location in the HR diagram 

In order to compare the observations with the theoreti- 
cal predictions, it is necessary to estimate the position of 
VV Ser in the HR diagram. As discussed in Section 2, spec- 
tral type and distance to VV Ser are quite uncertain. First, 
we have estimated the T e g value from the adopted spectral 



type A2±2 subclass by using the Schmidt-Kaler (1965) 



conversion tables, noting that for these relatively early- 
types these values are in agreement with the more recent 
estimates of, e.g., Kenyon & Hartmann (1995) (see their 



Table A5) . As a result, we obtain T eff =9000±1000 K. 
Concerning the luminosity, we decided to investigate all 
the values in the range 1.3 < logL/L Q < 2.3, covering 
the empirical estimates in the literature. The studied re- 
gion is represented by the large box in the HR diagram 
shown in Fig. El that includes a set of PMS evolution- 
ary tracks computed using the CESAM code IjMorel 19971 
|Marques et al. 2004| ) for stellar masses from 1.5 to 4.5 M Q , 
and the theoretical instability strip computed for the first 
three radial modes by Marconi fc Palla (1998)[ on the ba- 
sis of nonlinear convective pulsation models. We note that 
for the estimated effective temperatures of VV Ser we do 
not expect pulsation in the first three radial modes since 



the predicted instability strip is cooler at each luminosity 
level. We might instead expect pulsation in higher radial 
overtones and/or in nonradial modes. In fact the instabil- 
ity strip is predicted to move toward higher effective tem- 
peratures as the radial overtone increases (see Grigahcene 
et al. 2006). 

5.2. Oscillation frequencies 

In order to evaluate the sensitivity of the predicted pe- 
riodicities to the input stellar parameters, we computed 
a fine grid of structure models along the CESAM evolu- 
tionary tracks covering a mass range from 2.2 to 4.0 M Q 
and effective temperatures from ~ 8000 K to ~ 10,000 K. 
The physical properties of the selected PMS models are 
reported in Table El (available in electronic form). 

For each model, frequencies for 1=0,1,2 1 
modes were computed using the Aarhus adi- 
abatic nonradial pulsation code (ASTEC- 
http : / / astro . phys . au . dk/ ~ j cd/adipack . n/ ) and 
compared with the observed periodicities. The range 
covered by the predicted frequencies depends on both the 
stellar mass and the effective temperature, moving toward 
higher values as the stellar mass (effective temperature) 



1 The reason why we consider only 1=0,1,2 modes is that, as 
noticed by other authors (see e.g. Suran et al. 2001, Baglin et 
al. 2000), these modes are expected to be easily detectable, as 
their visibility coefficients remain large after integration over 
the whole stellar disk. 
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Fig. 5. PMS evolutionary tracks for the labelled masses 
and solar chemical composition (Z = 0.02, Y = 0.28) 
computed with the CESAM code. The dashed box rep- 
resents the estimated uncertainty on the position of VV 
Ser. It is also shown the location of the theoretical insta- 
bility strip for the first three radial modes predicted by 
Marconi & Palla (1998) on the basis of nonlinear convec- 
tive pulsation models. The empty circles represent the 32 
models run for the pulsational analysis. The filled symbols 
are for the best fit models discussed in Section 5. 



decreases (increases) at fixed effective temperatures 
(mass). In Figs. 16191 the predicted frequencies with 1=0 
(filled symbols connected by dotted lines), 1=1 (filled 
symbols connected by long dashed lines) and 1=2 (filled 
symbols connected by dashed lines) are compared with 
the empirical periodicities (labelled horizontal solid lines) 
for the labelled masses and effective temperatures. 

We find that the three lowest frequencies, namely /2, fz 
and /4, cannot be reproduced simultaneously by p modes 
in the explored mass and effective temperature range. In 
particular, / 2 is only consistent with a g mode (1=1,2 with 
n < 0) unless the effective temperature is much lower than 
the estimated empirical range. To illustrate this point, we 
show in the left bottom panel of Fig.[!2]the plot for a signif- 
icantly cooler model (T off =6997 K) with M=4.0 M . In 
this case, the range of the three lowest observed frequen- 
cies is covered with the predicted I = 0, n = 0, 1, 2 modes, 
with f-2 close to the theoretical result for the fundamental 
mode. Pulsation in these low order radial modes at such a 
low effective temperature is also expected on the basis of 
the location of the predicted instability strip (see Fig. |3J). 

A tentative mode identification based on the various 
plots is summarized in Table available in electronic 
form. Model frequencies are required to match the ob- 



served ones within 2.5/i-ffz, in order to take into account 
both the mean error on the measured periodicities (~ 0.5 
fiHz) and an estimate of the (unknown) model intrinsic 
uncertainty. The models reproducing all the observed fre- 
quencies with Z=0,l,2 p and/or g (or /) modes are the 
ones reported in Table 21 and shown in Fig. [S] as filled 
circles. The stellar mass and luminosity associated with 
these models range from 3.6 to 4.0 M Q and from 125.2 to 
207.6 Lq. To remove the degeneracy on the model solu- 
tion, an independent empirical estimate of the large sepa- 
ration would be very important. However, the evaluation 
of the large separation from the frequency power spectrum 
would require more accurate, long-time data. 



Table 4. Preliminary mode identification as a function 
of the explored model input parameters for models repro- 
ducing all the seven observed frequencies. The mass is in 
solar units (M Q ) and effective temperature is in K. 
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8 253 
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-3 
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We note that, among the best-fit models, mod37 is in- 
side the theoretical instability strip for the three lowest 
radial modes, so that fi and are correctly reproduced 
by pulsation in the fundamental and the second overtone, 
respectively. In this case all the seven observed frequen- 
cies can indeed be reproduced by p modes. Therefore, the 
results seem to suggest that either the effective tempera- 
ture found in the literature is significantly overestimated, 
or that we have to admit the possibility of pulsation in at 
least one g mode. 
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Available studies point towards the existence of unsta- 
ble low degree p-modes in PMS stars. Consequently, in 
the above analysis, models with p-modes in the observed 
frequency range have been preferred. The possibility that 
g-modes are also excited to the observed amplitudes has 
yet to be confirmed by a detailed theoretical analysis of the 
excitation mechanisms of such modes in young stars and 
supported by observational mode identification. A more 
detailed analysis (theoretical and observational) in this 
direction is clearly required. 



5.3. Effect of rotation on mode identification 

Here, we give a simplified, preliminary account of the ef- 
fect of rotation on nonradial modes using the asymptotic 
relation for the rotational displacement. This is given by 
mfrot, where m is the spherical harmonic number and 
frot = {vrot sin i)/(2ttR sin i) the rotational frequency. 
The term v ro tsmi can be taken from published values. 
According to Sect. 3, the rotational velocity varies be- 
tween a minimum of w rot sini ~85 km/s and a maxi- 
mum of 200 km/s. For the sini term in the denomina- 
tor, we need to make an assumption on the inclination 
angle i. For simplicity, we consider two extreme cases, 
namely i\ — 7r/4 and %2 = tt/2. Finally, for the theo- 
retical radius we take the values of our best-fit models 
as given in Table 01 Adopting v ro tS,uii =200 km/s, we 
find a rotational splitting 8.5(17) < A/i < 14.6(29.2)^Hz 
and 6(12) < A/ 2 < 10.3(20.6)^Hz for m = ±1(2), re- 
spectively. Conversely, for a velocity of 85 km/s the split- 
tings are reduced to 3.6(7.2) < A/i < 6.2(12.4)//Hz and 
2.5(5.0) < A/ 2 < 4.4(8.8)/zHz for m = ±1(2), respec- 
tively. We have also calculated the upper limit of Af under 
the extreme assumption of v < Vbreakup-, where Vbreakup is 
the break up velocity for the models reported in Table 
yielding Af rn ax ~ 20/iHz and thus suggesting that the 
rotational splitting might introduce significant error bars 
on the theoretical mode identification. 

On the other hand, if we assume that the measured 
periodicity at 2.2 d is associated with rotation, we derive 
a rotational splitting of ~ ±5(10)/ii/z for m = ±1(2), in 
better agreement with the values m x 2.5 to m x 4AfiHz 
obtained for a rotational velocity of 85 km/sec and a ra- 
dius range covered by the best fit models reported in Table 
0]and a star seen close to edge-on (i — tt/2). Thus, these 
initial results seem to indicate that rotation may affect sig- 
nificantly the theoretical mode identification, with larger 
radii providing smaller rotational splitting for a given ro- 
tational velocity. However, since rotation alters the struc- 
tural and evolutionary properties of stellar models, an ac- 
curate estimate of the effects on mode identification can 
only be obtained by a specifically designed numerical code 
that incorporates all the relevant physical ingredients (see 
e.g. Reese et al. 2006). 
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Fig. 6. Comparison between predicted (filled symbols) 
and observed (horizontal solid lines) periodicities for mod- 
els at 3.2 Mq (bottom panels), and 3.4 M© (upper panels). 
Symbols are connected by a dotted line in the case of /=0, 
by a long-dashed line for 1=1 and by a dashed line for 1=2. 
The two panels refer to the labelled effective temperatures. 
When a predicted frequency reproduces an observed one 
the corresponding filled symbol (along the given 1 slanting 
line) is located at the intersection with the horizontal line 
representing the observed value. 
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Fig. 7. The same as for Fig. but for 3.6 Mq (bottom 
panels) and 3.7 Mq (upper panels). 
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6. Conclusions 

We have detected S Scuti-type pulsation in the Herbig Ae 
star VV Ser on the basis of ground-based observations 
carried out during three consecutive years. In particular, 
using the best dataset obtained in 2004 seven pulsation 
frequencies have been measured with values ranging from 
2.69 to 10.24 c/d. These results allow to include VV Ser 
in the list of known PMS stars with 5 Scuti variations ( 



IZwintz et al. 20041 |Ripepi et al. 2006| ). In order to com- 
pare the extracted periodicities with model predictions, 
we have considered a large empirical range both in lu- 
minosity and effective temperature, reflecting the exist- 
ing uncertainties on the stellar properties. Within this 
space of physical parameters, we have computed a fine 
grid of inner structure models along the CESAM evolu- 
tionary tracks. The corresponding pulsation frequencies 
for /=0,1,2 modes have been evaluated by means of the 
ASTEC code and compared with the observed ones by 
requiring an agreement within 2.5 fj,Hz. 

As a result, we find a number of best fit models cor- 
responding to PMS stars with mass 3.6-4.0 M and lu- 
minosity log L/Lq w2.1-2.3 that reproduce the seven ob- 
served periodicities with radial and/or nonradial p and g 
modes. Such an independent confirmation of the funda- 
mental stellar parameters range previously estimated for 
VV Ser is an important byproduct of the main objective 
of asteroseismology (which will be addressed in a future 
work, requiring more and better observations and more 
accurate modelling), i.e. of sensing the internal structure 
of the star. We also notice that if the effective temperature 
is significantly lower (~ 7000 K) than the values based on 
the empirical spectral types (~ 9000 K), we are able to 
reproduce all the observed periodicities with p modes for 
a stellar mass of 4 Mq and a luminosity of 125 L Q . An ac- 
curate spectral type determination of VV Ser is urgently 
needed in order to discriminate among the different so- 
lutions. In the meantime, the present analysis seems to 
indicate a stellar mass in the range of the largest masses 
expected (and indeed observed) for this class of pulsators. 
The possible effect of rotation on the predicted frequen- 
cies depends on the adopted estimates of the rotational 
velocity and inclination angle, as well as on the stellar ra- 
dius. If the measured periodicity at 2.2 d is interpreted 
in terms of rotation, the corresponding rotational split- 
ting is of the order of m x 5/jHz. Finally, we note that 
further photometric observations of VV Ser, especially by 
means of multisite campaign, would be of great help to 
better define the frequency spectrum and to eliminate the 
uncertainty on the observed frequencies due to the 1 c/d 
alias. 
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Table .1. Log of the observations. 



HJD-2450000 

start (days) 
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end (days) 


Duration 
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6.1 


53185.350 


53185.597 
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53186.420 
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53209.544 
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53210.399 


53210.518 


2.9 



Table .2. Physical properties of the selected PMS models. 
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202.2 


9376 


1.22 


mod32 


3.8 


4.97 


231.2 


10106 


1.26 


mod33 


3.9 


6.39 


166.5 


8210 


1.06 


mod34 


3.9 


6.00 


190.7 


8766 


1.09 


mod35 


3.9 


5.59 


218.1 


9392 


1.13 


mod36 


3.9 


5.16 


248.4 


10092 


1.16 


mod37 


4.0 


7.63 


125.2 


6997 


0.89 


mod38 


4.0 


6.58 


180.5 


8253 


0.98 


mod39 


4.0 


6.20 


205.2 


8760 


1.02 


mod40 


4.0 


5.80 


232.9 


9366 


1.05 


mod41 


4.0 


5.40 


263.8 


10023 


1.08 
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Table .3. Preliminary mode identification as a function of the explored model input parameters. The mass is in solar units 
(Mq) and the effective temperature in K. 



Model 

Mass 


Toff 


I 


/i 


.ft h .A h 

Possible value of n 


h 


h 


modi 

2.2 


8 749 



1 

2 




-5 


-3 






mod2 

2.2 


9 913 



1 

2 


-5 


-7 


-4 -5 
-7 -9 


-2 
-4 




mod3 

2.4 


9 304 



1 

2 




-5 
-9 


-3 

-5 -6 






mod4 

2.4 


10 000 



1 

2 




-6,-7 


-4 
-8 




-2 


mod5 

2.6 


8 459 



1 

2 






1 


1 

-1 
-2 




mod6 

2.6 


9319 



1 

2 


-1 


mod7 

2.6 


10023 



1 




mod8 

2.8 


8168 



1 


1 

-1 
-2 




3 




2 
2 


mod9 

2.8 


9 885 



1 




modlO 

3.0 


8 007 



1 

2 


-1 


-4 
-7 


-2 

-4 3 




4 


modll 

3.0 


8 758 



1 

2 


1 

-1 
-2 


-5 
-8,-9 


3 

-2 

-5 




2 
2 


modl2 

3.0 


9 621 



1 

2 


-3 


-6 
-10 


-5 -6 






modl3 

3.0 


10 547 



1 

2 


-4 


-6 
-11 


-3 -4 
-6 -7 






modl4 

3.2 


8 573 



1 

2 


-1 


-4 
-7 


-2 

-3 -4 3 




4 


modl5 

3.2 


9417 



1 

2 


1 

-1 
2 


-5 
-9 
-8 


3 

-2 

-4 -5 




2 
2 


modl6 

3.2 


10 375 



1 

2 


-3 


-5 
-10 


1 

-5 -6 


1 




modl7 

3.4 


8 086 



1 

2 


1 


-3 
-5 


1 5 
-1 

-2 -3 


2 




modl8 
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Table .3. continued. 



IVloael 






£ 

h 


h 


h 


h 


£ 

h 


£ 

h 


77 


Mass 


rr\ 

J-cS 


7 
( 




Possible value of 


n 






3.4 


8 762 


1 




-3 






4 




4 






A 




-6,-7 


-3 


-4 








modl9 

















2 




3.4 


9 538 


1 




-4 








1 








A 




-7,-8 












mod20 





















3.4 


10414 


1 




-5 




-3 












A 




-9 






1 


-1 




mod21 





















3.6 


8 234 


1 








-1 


5 










A 


1 
i 


-5 


-1 


-2 




Q 
O 




mod22 







2 










3 


5 


3.6 


8 868 


1 


1 


-3 


-1 














A 




-6 


-2 


-3 




o 
z 


4 


mod23 















4 




4 


3.6 


9 590 


1 








-2 












A 




-7 


-3 


-4 


■j 
o 


1 
1 




mod24 







1 














3.6 


10402 


1 




-4 


-2 








2 






A 


-A 


-8 


-4 


-5 


o 

z 


U 




mod25 







3 






1 


6 


4 


6 


3.7 


8 211 


1 




-2 




-1 












A 


o 
A 


-4 


-1 


-2 


F 



Q 

o 




mod26 











1 




5 






3.7 


8 823 


1 


1 


-3 


-1 






2 


5 






A 




-5 


-2 


-3 








mod27 





















3.7 


9 518 


1 




-3 






4 










A 


(J 


-6 












mod28 

















2 




3.7 


10 300 


1 




-4 




-2 




1 








A 


-1 


-7,-8 




-4 








mod29 







3 






1 








3.8 


8 155 


1 


2 


-2 






6 


4 


6 






A 


o 

A 


-4 





-1 








mod30 





















3.8 


8 727 


1 








-1 


5 










A 


1 


-5 


-1 


-2 




Q 

o 




modal 







2 


-3 








3 


5 


3.8 


9 376 


1 


1 


-6 


-1 














A 




-5 


-2 


-3 




o 
z 


4 


mod32 















4 




4 


3.8 


10 106 


1 








-2 












z 




-7 


-3 


-4 


Q 
O 


1 




mod33 











2 










3.9 


8210 


1 


2 


-2 






6 


4 








2 




-4 





-1 






6 


mod34 


















4 


6 


3.9 


8 766 


1 




-2 




-1 












2 




-4 


-1 


-2 


5 


3 


5 


mod35 











1 




5 






3.9 


9 392 


1 


1 


-3 


-1 






2 








2 




-5 


-2 


-3 


4 
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Table .3. continued. 



Model 






h 


h 


h 


/4 


h 


h 




1\ Jnnn 

Mass 


Jeff 


7 

I 




Possible value of 


n 






mod36 





















3.9 


10092 


1 

A 


U 


-3 
-6 










4 


JOT 

mod37 







5 


1 


3 




9 




9 


4.0 


6 997 


1 








1 




6 








A 


4 


-2 


2 


1 


o 




8 


mod38 











A 








7 


4.0 


8 253 


1 


2 


















2 




-3 





-1 


6 


4 


6 


mod39 







3 






1 


6 


4 




4.0 


8 760 


1 

2 


2 


-2 
-4 




-1 






6 


mod40 











1 










4.0 


9 366 


1 

2 


1 


-5 


-1 


-1 

-2 


5 




5 


mod41 







2 










3 




4.0 


10 023 


1 




-3 


-1 














2 





-6 


-2 


-3 




2 


4 



